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Abstract 

The article deals with a number of the existing variants of direct calculation 
of amplitudes of processes with polarized Dirac particles. It is shown, that all of 
them are special cases of one and the same mathematical scheme. The advantages 
and disadvantages of this scheme are considered. A new variant of the method of 
calculation of amplitudes, keeping all the advantages of this scheme, but free from 
disadvantages, is proposed. In particular, this variant is suitable for the evalua- 
tion of amplitudes of processes with interfering diagrams and does not need any 
additional calculations, that are necessary in the general case. Expressions for the 
amplitudes of processes involving both massive and massless particles are presented. 
These expressions make it possible to create in an easy way the computer programs 
for automatic calculation of amplitudes. The existing computer programs for cal- 
culation of amplitudes, their limitations and disadvantages are briefly considered. 



1 Introduction 

In order to calculate the observables of the processes with the Dirac particles, especially 
for high-order diagrams and in the case of taking into account the polarization effects, it 
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is necessary to calculate the traces of products of great number of the Dirac 7-matrices. 
Thus, we face the problem of obtaining the analytical expressions for the physical quan- 
tities we are interested in. 

One way to eliminate this problem is to calculate the amplitudes of the processes di- 
rectly. This method allows as to reduce essentially the number of the Dirac 7-matrices 
in the expressions considered so also the number of the expressions themselves (in the 
situation when we have the interfering diagrams). The results obtained in this way de- 
pend only on 4-vectors and are suitable for numerical analysis. To calculate the squared 
amplitude (by means of which observables are expressed), we need to calculate only the 
squares of the modules of the complex numbers obtained. The numerical values obtained 
in this way are equivalent in principle to the values obtained in the classical way, but can 
have smaller errors caused by loss of accuracy at the calculations. 

The appeal of such an approach has stimulated a lot of articles, where a number of 
different ways of calculation of amplitudes have been proposed. However all of them have 
the limitations. Moreover, there are the incorrect methods (for example [[FjJ [35], see 
subsection 6 of Section 6). These reasons and the fact that the expressions which have 
been obtained by different ways for the same amplitudes have different analytical form 
lead to a certain mistrust to the idea of direct calculation of amplitudes. 

Up to now the methods of calculation of amplitudes have been separated into two 
classes: the covariant methods, i.e. the methods reduced to the evaluation of the traces 
of the products of the Dirac 7-matrices (see jl| [p5|D and the methods reduced to the 
multiplication of 7-matrices and bispinors which are written in the matrix form (see [S3 - 
|36| , |38[1 f42fl). However all of them are in fact the special cases of one and the same 



mathematical scheme. 

In Section 2 the covariant methods of calculation are briefly considered. The general 
mathematical scheme of calculation of amplitudes is presented. It describes all the known 
methods. It is shown also how the results of any special method can be transformed into 
the results of another one. 

In Section 3 the advantages and disadvantages of the general scheme are considered. 
In particular we consider carefully the additional calculations which are necessary to be 
performed in the general case when one has the interfering diagrams. 

In Section 4 a new covariant method of calculation of amplitudes is proposed. This 
method keeps all the advantages of the general scheme but is free from its disadvantages. 

In Section 5 the general expressions for the amplitudes of processes involving both 
massive and massless Dirac particles are presented. (Examples of calculation of the con- 
crete processes can be found in 0.) Using the presented formulae one can create the 
computer program for automatic calculation of amplitudes. 

In Section 6 the methods of calculation of amplitudes based on the multiplication of 
7-matrices and bispinors which are written in the matrix form and some similar methods 
are considered. It is shown that all of them are the special cases of the general scheme 
considered. 
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In Section 7 the existing computer programs for automatic calculation of amplitudes, 
their limitations and drawbacks are briefly considered. 



2 The general scheme of the covariant calculation of 
amplitudes 

There is an even number (2N) of fermions in initial and final state for any reaction 
with Dirac particles. Therefore every diagram contains iV nonclosed fermion lines. The 
expression 

M if = u f Q Ul (1) 

corresponds to every line in the amplitude of the process, Uj, Uf are the Dirac bispinors 
for free particles. (For definiteness, we assume that fermions are particles. However the 
results obtained can be generalized to the case when both of the fermions are antiparticles 
or one fermion is a particle and another one is an antiparticle. This generalization is 
considered at the end of Section 4.) 

u = n + 7° , 

and Q is a matrix operator characterizing the interaction. The operator Q is expressed as 
a linear combination of the products of the Dirac 7-matrices (or their contractions with 
4- vectors) and can have any number of free Lorentz indexes. 
The formula ([[]) can be rewritten as 

M if = Tr(QuiU f ) 

however the expression for the operator UiUf is not known. Because of this, in order to 
calculate Mif we use the following general scheme 



M if = UfQui = (ufQui 



UiZuf Tr(QuiUiZufUf) 



UiZuf UiZuf 
Tr(QuiUiZufUf) _ Tr(QuiUiZufUf) 



(2) 



if 



\ u i^ u f\ ^Tr(ZuiUiZufUf) 
where Z is an arbitrary 4 x 4-matrix, 

Z = 7°Z + 7° 

(symbol ~ stands for "an equality up to a phase factor sign"). 

The projection operators are substituted for uu in (|]). For a particle with mass m we 
have: 

u(p, n)u(p, n) = -(p + m)(l + 75^) = V{p, n) (3) 



whereQ a = 7 M a M for any 4-vector a, 

n i n o 

p = m , n = —1, pn = , uu = 2m , 75 = i'y 7 7 7 

(p is the 4-momentum of the particle, n is the 4-vector specifying the axis of the spin 
projections of the particle). 

For a massless particle the projection operator is the following 

u±{p)u±(p) = -(1 ± 7 5 )p = V±(p) (4) 

where 

P 2 = , m±7a 1 «± = 2p M 

(signs ± correspond to the helicity of particle) . 

Actually the general scheme (0) describes all the known methods of calculation of 
amplitudes. 

1. In the articles [|J - |§ one chooses 

Z= 1 . 



The results of the articles [0 reduce to the same choice for the calculation of 
amplitudes of type u±(pf)u T (pi) for the processes with massless particles. 

2. In 0, M one proposes 

Z = 75 • 



The results of the article 11 reduce to the same choice. 



3. The results of Ul2| , [ |T3p reduce to 

Z = l° • 

The expressions obtained in the article [frjj reduce to the same choice under the 
following restrictions: 

(a) Dirac particles with equal mass in the center-of-mass frame are considered; 

(b) the polarization state of the particles is the helicity, i.e. 4-vectors, which deter- 
mine the axes of spin projections, are (see e.g. |60| ): 



1 We use the same metric as in the book ]3j : 

a M = (ao, a), a M = (ao, —a), ab = a^U L = a^b® — ab 
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(c) the particles have the opposite helicities. 
4. The results of |L4| also reduce to the choice 

Z = 7°7 5 Z 

1 



[where ¥ 



0, PxPz, PyPz,~(Px +P y ) ,1 = -1, Oi) = (/p/) = 0] 



\P\\JPx 2 +Py 2 

under the following restrictions: 

(a) Dirac particles with equal mass in the center-of-mass frame are considered; 

(b) the polarization state of the particles is the helicity; 

(c) the particles have equal helicities. 

5.1. In [|l|] it was proposed to choose for Z 

Z = I 

where I is an arbitrary 4- vector such that I 2 = —1. The same choice was proposed 
in H for the calculation of amplitudes of processes with massless particles. 



5.2. In [TB] it was proposed to choose for Z 



(where q is an arbitrary 4- vector such that q 2 = 0, it was proposed to choose 
q^ = (1,1,0,0) for numerical calculations) under the following restrictions: 

(a) the 4-vectors, which determine the axes of the spin projections, are 

1 



n = A 



m 

m ypq) 



A = ±1 



(p is the 4-momentum of the corresponding particle), 



(b) 



A,; = A 



/ 



5.3. In 0| it was also proposed that for the calculation of amplitudes of processes with 
massless particles 

Z = r 



where r is an arbitrary 4-momentum such that r 2 



m 
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6. The results of papers JT7J] reduce to 

Z=l+ 7 ° . 



The results of the paper |20[ reduce to the same choice under the restriction 

0/)o 



n f 



v n/) , ft* - — — — p. 



( n f)o 



(p-do + m (p f ) + m f 



(Pffk) ~ ( J 1 *® — (PfPi 



Pf 



(pi)o + mi 



7. The results of []20| also reduce to 



where 



Z= (l + 7°)7 5 / 
-(n/) , -n* + 



(ft;Jo 



P 



i n f\ 



m f 



(pi)o + 0/) + m/ 



P/ 



(p»)o + m i 



where 



F = (0,-o), o 2 = l, {ofii) 



(ft;Jo 



(Pi)o + m. 



(op*) • 



8. The results of papers [^T|], |22| reduce to 

Z = m + f 

(r is an arbitrary 4-momentum such that r 2 = m 2 ), as this takes place in this paper 
for the 4-vectors, which determine the axes of the spin projections, one uses 



ft; 



m hi - (PiPf)Pi 



{PiPff 



mfrrif 



n .f 



m }Pi ~ (PiPf)Pf 
m f\J{PiPf) 2 ~ m^nif 



9. The results of |23j for the amplitudes of type u±(pf)u T (pi) with massless particles 
reduce to 

Z = qiq 2 

where q±, q 2 are arbitrary 4-vectors such that q\ = q\ = 0. 
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10. The results of the papers |2I|| , |25| ] for the amplitudes of type u±(pf)u T (pi) with 
massless particles reduce to 

Z = lq 

where I, q are arbitrary 4- vectors such that 

P = -l, q 2 = 0, (ql) = 0- 



As this takes place one chooses P = (0,0,1,0), g M 
calculations. 



^1,1,0,0) for numerical 



The results of the paper |16| reduce to the same choice for the calculation of ampli- 
tudes of processes with massive particles under the following restrictions: 

(a) for the 4-vectors, which determine the axes of the spin projections, one uses 

1 



n = A 



m 

—P - 

m ypq) 



A = ±1 



(b) 



(p is the 4-momentum of the corresponding particle), 



A,- — — A/ 



The remaining methods of the calculation of amplitudes are considered in Section 6. We 
emphasize again that all this methods can be reduced to (0). 

The results of the calculation of amplitudes for different Z differ only by the phase 
factor. Really, 



M if {Z x ) 



TrjQViZiVf) TrjQVjZrVf) Tr(Z 2 ViZ 2 V f ) 
^Tr(Z{P t Z{P f ) ^Tr(Z{PiZxP f ) Tr{Z 2 ViZ 2 V f ) 

Tr{QViZ{P f ) ■ Tr{Z 2 ViZ 2 V f ) 



Tr{QV l Z 2 V f ) ■ Tr(Z 2 V i Z 1 V f ) 



(5) 



Tr{Z x ViZ 2 V f ) ■ Tr(Z 2 ViZ{P^Tr{Z 2 ViZ 2 Vi) 



Tr{Z 2 ViZ 2 V f 



TrjZiPiZiV} 
Tr^ViZ^j 



1/2 



M if (Z 2 ) ■ e 



i(f> 



where 



J4> 



TrjZzViZ^f) 
Tr{Z{PiZ 2 V s ) 



1/2 



(6) 
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In (|5|) we have used the identity: 

Tr{AViBV f ) ■ Tr(CViDV f ) = Tr{A Ui UiBu f u f ) ■ Tr(CuiUiDu f u f ) 

= (u f Aui)(uiBuf)(ufCui)(uiDu f ) = (u f Aui)(uiDuf)(ufCui)(uiBu f ) (7) 

= TriAuiUiDujUf) ■ Tr(CuiUiBufUf) = Tr(AViDV f ) ■ Tr(CViBV f ) 
where A, B, C, D are arbitrary 4 x 4-matrices. 

3 Advantages and disadvantages of the general scheme 

Let us consider the advantages and disadvantages of the general scheme of calculation of 
amplitudes (Q) illustrating them by the elementary example with Z = 1 . 
The advantages are following: 

1. The calculation of the amplitude of the process is analogous to the calculation of the 
squared amplitude of this process 

\ M if\ 2 — (u f Qui)(u f Qui)* = (u f Qui)(uiQuf) = Tr(QuiUiQu f Uf) = Tr(QViQVf) 

(8) 

and reduces to the standard operation of the evaluation of the trace of a linear 
combination of the products of the Dirac 7-matrices. 

As this takes place, for the calculation of amplitudes it is enough to make the 
replacement in (|8|) 

Q — Z (9) 

^Tr(ZViZVf) 

This observation allows us to supplement the computer algebra systems, using stan- 



dard methods for calculation of cross sections (e.g. fS3fl), by the programs for 
the calculation of the amplitudes. 

For Z = 1 the replacement (^j takes the form 

Q - 



yj[mim f + (piPf)} [1 - (ninf)] + ipin^ipfUi) 

2. The expression for the amplitude of the process contains lesser number of the Dirac 
7-matrices than the expression for the square of the modulus of the amplitude of the 
same process. Because of this, after the evaluation of traces of the 7-matrices the 
expression for the amplitude contains lesser number of terms than the expression for 
the square of the modulus of the same amplitude. The estimate of a gain for Z = 1 
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is presented in Table [TJ. (Here we consider the processes involving massive particles. 
The maximum number of the 7-matrices in projection operators of the massive 
particles is contained in the terms pa^hi and . Therefore the maximum 

contribution from the projection operators in expressions considered is equal to 4 
7-matrices.) 



Number of 7-matrices in Q , N 


1 


3 


5 


Number of 7-matrices in the term of M^/ 
maximum contribution, N± = 4 + 2N 


2 , which gives the 


6 


10 


14 


Number of the terms after the trace evaluation, (Ni — 1)!! 


15 


945 


135135 


Number of 7-matrices in the term of M.if , which gives the 
maximum contribution 2 , N 2 = (4 + N) — 1 


4 


6 


8 


Number of the terms after the trace evaluation, (N2 — 1)!! 


3 


15 


105 


(Ni - 1)!! 
The estimate of a gain, — 


5 


63 


1287 



Table 1: The estimate of a gain obtained at calculation of the amplitudes for Z = 1 

It is follows from Table [I] that the more complicated the process considered (i.e. the 
greater number of the 7-matrices is contained in the operator Q which describes 
the interaction) the greater gain is given by the method of the calculation of the 
amplitude in comparison with the standard method of the calculation of the square 
of the modulus of the amplitude. This is also true for any operator Z in the general 
scheme (Q). 

Really, when the number of 7-matrices in operator Q increases by /, their number 
in the numerator of the expression (Q) for the amplitude increases only by /, (de- 
nominator does not change), but in the expression (||) for the squared amplitude the 
number of 7-matrices increases by 21 (as the formula (§) contains both the operator 
Q and the operator Q). 

However, the scheme considered has two essential disadvantages in the general case: 

1. There is a denominator in (^), hence the ambiguity of the type - can appear during 

the calculations. For example if one chooses Z = 1 then the expressions for the 
amplitudes have the denominator 

y/Tr[P(pi,ni)V(pf,n f )] = \f\mmf + (piPf)} [1 - JjhJif)} + (pin f )(p f ni) 
which is equal to zero for 

n f = -rii + — -{pi + — p f ) 

rriimf + [PiPf) m f 



2 In this case (4 + N) is an odd number, and consequently the term with (4 + N) 7-matrices vanishes 
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2. All the expressions for the amplitudes [as it follows from (Q)] are known up to a phase 
factor: 

Mif = M if ■ . (10) 



\u 



■Zu f \ 



Obviously this fact creates no problem, if we calculate the amplitude for alone 
diagram. Really 



(M if y 



[Tr(QViZV f )Y [TriQmUiZufUf)]* {{ufQu^UiZuf)]* 



Tr(ZViZVf) sjTr(ZViZVf) pr(ZV % ZV s ) 
UfZui)(iiiQuf) _ Tr(ZuiUiQufUf) Tr(ZViQVf) 



111 



Tr(ZViZVf) 

M if (M if 



Tr(ZViZVf) sjTr(ZViZV f ) 
TriQVtZVf) ■ Tr{ZViQVf) 



TriQV&Vf) ■ TriZVtZVj 



Tr(ZViZVf) 

Tr(QViQVf) = |M, 



(12) 



Tr(ZViZVf) 
[In (0) we used the identity (0).] 

However, if we have the interfering diagrams, then the application of (|2|) leads to 
the fact that the expressions for the amplitudes corresponding to different channels 
are multiplied by different phase factors in the general case. In this situation the 
additional calculations are necessary. These calculations are considered below. 

Let us consider the diagrams of the process of the general form, which proceeds in two 
different channels (see Fig. [l]). 




Figure 1: The general form of the diagrams of the process proceeding in two interfering 
channels 



The expression 

M = (ufQm) ■ (ufRui,) = M if ■ M vr , 
corresponds to the first diagram. The expression 

W = (iifSui) ■ (u f Tu v ) = M ir ■ M vs 



(13) 
(14) 
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corresponds to the second one, where Q, R, S, T are arbitrary matrix operators charac- 
terizing the interaction. 

There are three possibilities to calculate the amplitude of this process correctly: 

(a) One can apply the Fierz transformation for the second diagram: / <-> /'. However 

this method requires a large number of additional calculations. 

(b) One can multiply both of the amplitudes by the same phase factors, for example by 

UiZiUf UiiZ 2 Ufi 
\uiZ\Uf\ \ui>Z 2 Uf>\ 

In this case we have the expression for the first diagram 

TrjQViZiVf) Tr(RV v Z 2 V r ) 
y/T^zTP&Vf) ^Tr[Z 2 V v Z 2 V r ) 

and for the second diagram 

Tr(SV l Z 1 V f TV il Z 2 V f ) 
yjTriZiV&Vf) ■ ^Tr(Z 2 V v Z 2 V r ) ' 

It is obvious that the calculation of the amplitude for the second diagram is com- 
plicated enough. Besides, it is not convenient to use the expressions which have 
different structure for the calculation of amplitudes for the different diagrams. 

(c) The third possibility is to calculate the relative phase for the first and second diagrams 

and to use the expression obtained for the phase correction of one of two amplitudes. 

We will use the third possibility. Interference term has the form 

M ■ (M'Y = M if ■ Mvp ■ (M if ,)* ■ (M ilf )* = {ufQuiXufRueXuiSuf^fuf) . 
Let us multiply the interference term by 

(u i Z 1 Uf)(u il Z 2 Uf>)(ufZ 3 u i )(ufZ 4 Ui>) (ufZ 1 u i )(uiZ 3 Uf)(ufZ 2 u i i)(ui>Z 4 Uf) = 

{u i Z 1 Uf){u i 'Z2Ufi)(Uf'Z 3 Ui){UfZ 4 Uii) (UfZ 1 U i )(UiZ 3 Uf')(Uf'Z 2 Ui')(Ui'Z4Uf) 

where Zi, Z 2 , Z 3 , Z 4 are still arbitrary matrix operators. 
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As a result in this case we have 
M-(M'Y 



Tr(qPiZiP f ) Tr{KP i >Z 2 V f >) 



Tr{Z{PiZ x Vf) Tr(Z 2 P v Z 2 P r ) 

X Tr(Z 3 pZ 3 P r ) ' Tr(Z,P t ,Z,P f ) ' T^pZ^Z^Z^) 
Tr{QViZ{P f ) Tr(RP v Z 2 P r ) Tr{Z 3 PSP r ) Tr{Z A P v TP f ) 



(15) 



TrWiZ&j) sjTr(Z 2 VirZ 2 V f ) y/Tr(Z 3 PiZ 3 P r ) ^Tr(Z 4 P if Z 4 P f ) 

Tr(Z 1 P i Z 3 P r Z 2 P ll Z 4 P f ) 



x 



Tr{Z 1 V i Z 1 Vf)Tr(Z 3 ViZ 3 Vf>)Tr(Z 2 V i >Z 2 Vf>)Tr(Z 4 V v Z 4 P f ) 

= M if ■ Mi>f> ■ (Mif)* ■ (Mi'/)* ■ K 

where M-if, M-vp, (MifY, {M-i'f)* are given by expressions analogous to (0) and (JTTj) , 
and the coefficient K is given by 

K = Tr(Z 1 P i Z 3 P fl Z 2 P i ,ZiP f ) 

/Tr(Z 1 P 4 Z 1 P / )Tr(Z 3 P l Z 3 P f )Tr(Z 2 ^Z 2 P / 0Tr(Z4^Z4P / )' 



It is obvious that 

\K\ = I. 

Thus we have to calculate the amplitude of process with interfering diagrams as follows 

M + M' = K-MifMi'f + MifMi' f . (17) 

In particularly, if 



Z\ — Z 2 — Z 3 — Z 4 — 1 

we have 

K = 



Tr{ViPf>Vi>P f ) 



lTr{P i P f )Tr{P i P r )Tr{P il P r )Tr{P il P f ) 
In the massless limit we have 

u ± (p f )Qu T (p^ TrmMlTl5)] 



K _ Tr]pip r pi>p f [l =F 7s)] ^ 
8 yj (PiPf) (piPf')(Pi>Pf) (Pi'Pf) 

The formulae ([jj§), C PL9|) generalize the method of the calculation of amplitudes proposed 
in PI. 
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4 A new method of the calculation of amplitudes 

All the disadvantages of the general scheme of the calculation of the amplitudes listed in 
the previous section can be eliminated with the help of a successful choice of the operator 
Z. 

As a suitable choice we propose 

Z = V [see (§)] or Z = V± [see (§] . 
Really, in this case the unknown phase factor for each line of the diagram splits into two 
parts [see (|10|)1: 

UiVUf , , UiUUUf , , UiU UUr 



\u{PUf\ \UiUUU,\ \U{U\ \uu,\ 

Let us consider again the diagrams of the process shown in Fig. |I[ 
Calculating the amplitudes (0), (0J), we have 

M = M if ■ M it f , ■ 
M' = M if i ■ Mi,, ■ 



UiU 


UU, 


Ui'U 


UU,' 


\u(a\ 


\uu,\ 


\u~i'U\ 


\uu,/\ 


UiU 


UUf, 


UiiU 


UU, 


\UiU\ 


\uu,i\ 


\Ui'U\ 


\uu,\ 



that is the amplitudes of the different diagrams are multiplied by the same phase factor, 
and therefore one can ignore it. This conclusion is right for any number of the interfering 
diagrams and for any number of the nonclosed fermion lines in the diagrams. 
We also note that, as would be expected [see (|TG|)1 

K = 1 

if 

Z\ = %2 — z% = z^ = *p 

or 

Z\ = ^2 = z% = z A = v± 

since projection operators have the following properties: 

V = V , VAV = Tr\VA] ■ V , (20) 
P± = V± , V±AV± = Tr[V±A] ■ V± . (21) 

Really 

p = 7°P + 7° = 7°(^) + 7° = 7 (m/ + 7°) + 7 

= 7°[(7°) + (« + ) + m + ]7 = 7 [7 Mn + ]7 = mi + 7° = uu = V , 
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TAP = {u) a {u)p{Afp{u) p {u) s = [{u) p {A)to(u) p ](u) a {u) s 

= [{u) p {u)p{Af%u) a {u) 5 = Tr[VA] ■ V . 

Therefore we can calculate the amplitude of the process with interfering diagrams as 
follows 

M + M' 



Tr(QV t VV f ) 


Tr(RV v VV r ) 


^Tr(VVi)Tr(VV f ) 


sjTr{VVi)Tr{VV r ) 


Tr{SV{PV r ) 


Tr(TVi'VV f ) 



lTr{VVi)Tr{VVf) ^Tr{VV v ) Tr (VV f ) 
Tr(QViVVf) ■ Tr{RVi<VV r ) + Tr{SV{PV r ) ■ Tr{TV v VV f ) 



[22] 



lTr{VVi)Tr{VVi>)Tr{VVf)Tr(Wfi 
[where for the denominators in we used the identity 



Tr(VViVV f ) = Tr{VVi)Tr{VV f ) , 

which follows from © - ©]. 

This expression enables us to calculate the amplitude numerically. The complex num- 
bers obtained in this way may be used to calculate the squared amplitude. 

Note that in and in the analogous formulae used further we shall use the equality 
sign instead of the symbol ~, since now there does not exist any trouble with the phase 
factors. (However one can not do this when the phase factor is not the same for all terms! 
- see, for example, subsection 5 of Section 6.) 

Thus, if we use the projection operator for Z then in addition to the two advantages 
of the general scheme (|2|) considered in Section 3 we obtain one more advantage: 

3. For the processes with interfering diagrams the number of the expressions which are 
necessary to be calculated is reduced since it is not necessary to calculate the inter- 
fering terms. 

If in individual cases under numerical calculations the denominator in (B^) vanishes, 
it is sufficient to change an arbitrary projection operator V entering ( ^2[) (i.e. to change 
arbitrary 4- vectors entering it). In this case the amplitudes for the different diagrams are 
multiplied by one and the same phase factor. Really, under the replacement of V by 
V' , we obtain [using (H) and 



TriQV.VVf) Tr(QViV'V f ) 

JVlif - 



Tr(VVi)Tr(VV f ) JTr(V'Vi)Tr(V'V 



Tr{V'ViVV 



fj 



Tr(VViV'V f ) 



1/2 
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Mi'f 



Tr{RVeVVf 



Tr(RVi>V'Vi 



f 



i'f 



Tr(VV % >)Tr(VV r 

Tr{SV{PV f ) 
Tr(VVi)Tr(VV f > 

Tr(TV v VV f ) 



lTr{VV v )Tr{VVf) 
As this takes place we have [see (0)] 

Tr(V'ViVVf) Tr(V'Vi'VVf) 



Tr(V'Vi>)Tr(V'Vf) 

Tr{SV{P'V fl ) 
Tr(V'Vi)Tr(V'V r ) 

TriTVyV'Vf) 
Tr(V'Vi>)Tr(V'V f ) 



Tr{V'VvVVf 
Tr(TVi>V'V f 



Tr{V'ViVV 



f'j 



Tr(VViV'V f ) 

Tr(V'Vi.W f ) 
Tr(VVi>V'V f ) 



1/2 



-,1/2 



1/2 



_Tr(VViV'V f ) Tr(VVi>V'V 



1/2 



Tr{V'ViVV f >) Tr{V'Vi>VV f 



Tr(VViV'V f >) Tr{VV v V'V 



1/2 



Note that if the operators, characterizing the interaction, contain the product of 7- 
matrices the number of which is greater than the number of 7-matrices involved in the 
projection operator, then the calculation of the amplitude for a single diagram is easer 
than the calculation of the squared amplitude. 

However, for the processes with interfering diagrams the method of the calculation of 
amplitudes is easier in any case, because we need not calculate the interference terms. 

As it was already mentioned, the method considered can be generalized easily to 
the case of the reaction with participation of antiparticles. To do this it is sufficient to 
substitute the projection operators of antiparticles in place of the operators of particles 
in (fj). For example, if we are interested in the value VfQui , where Vf is the bispinor of 
a free antiparticle, then 

Tr(QuiUiZvfVf) 



where 



VfQui 



v(p, n)v(p, n) 



p = m 

for a massive antiparticle, or 



n 



Tr(ZuiUiZvfVf) 

-(-m + p)(l + 7 5 n) 
, pn = , vv = - 



(23) 



(24) 



-2m 



v±(p)v±(p) = -(1=F7b)p 



(25) 



V 



, v±7/^± = 2 ^ 



for a massless antiparticle. 

As always, we use (^) or (f|) instead of Z. 
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5 The formulae for the calculation of the amplitudes 
of the processes involving the polarized Dirac par- 
ticles 

In this Section the expressions for the amplitudes of the processes with both the massive 
and massless Dirac particles are presented. Using the formulae (0) - (f|) and (^) - fl25|), 
we obtain: 

u±(p f )Qu±(pi) = u±(p f )Qv T (pi) = v T (p f )Qu±(pi) = v T (p f )Qv T (pi) 
_ Tr[Qpiqp f (l =p 7s)] ( 26 ) 

Here 

Z = ^(lTl 5 )q = V Tl q 2 = 0. 

The massless 4- vector q can be arbitrary, but it must be the same for all nonclosed fermion 
lines of the diagrams considered. 



u±(p f )Qu T (pi) = u±(p f )Qv±(pi) = v T (p f )Qu T (pi) = v T (p f )Qv±(pi) 
Tr[Qpi(m =p rl)pf(l =F 75)] 



(27) 



Here 



4y / [(rp i ) ± m{lpi)}[{rpf) =F m{lp f )} 
1 



Z = -(m + f)(l + 7 5 Z) =V, r 2 = m 2 , I 2 = -1, rl = . 



The same remark as for the vector q in (^) is takes place for the 4-vectors r and /. 

In the latter case we can not use the easier operator V± for Z, since in this case the 
numerator and the denominator are equal to zero. However to simplify the calculations 
we can require m = 0, r 2 = in (f27|): 

u±(p f )Qu T (pi) = u±(p f )Qv±(pi) = v T (p f )Qu T (pi) = v T (p f )Qv±(pi) 

7Y[Q&r%(l=F75)] (28) 
4:yJ(rpi)(rp f ) 

This formula generalizes the method of the calculation of amplitudes proposed in f2"l 



2q| . Further 

- ( \n ( \ - ( \n ( \ Tr [Q( m i ± PA + pi ±m i n i )qp f (1 =F 7s)] 



4y2[(gpi)±mi(gni)] (qp f ) 

(29) 
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- , \n ( \ - r \n i \ Tr lQ(- m i±Pini+PiTrn i h i )qp f (lTj5)} 
u±{pf)Qv{pi,ni) = v T {p f )Qv{pi,ni) = 



4 V 2 Km) T m(qni)} (qp f ) 

(30) 

«(p„ n f )Qu ±{pi ) = fi(p/ , n,)^) = Tr[9(1± /f M 7 T /^+^^^ )] , 

(31) 

-/ ^ / x _/ w / x Tr [Q(l ± j5)piq{-m f T Pfhf + Pf T m f n f )} 
v(p f , n f )Qu ± { Pi ) = v{p f , n f )Qv T { Pi ) = ===== , 

^2{qPi) [(qpf) =F mf(qnf)} 

(32) 

u(p f ,n f )Qu(pi,ni) 

_ TV [Qirrij ±pjJH +pj± mjhi)(l =F lb)qijn f Tp/Uf +Pf± m f h f )] ( 33 ) 
%\J[(QPi) ± ± rrif{qn f )} 

u{p f ,nf)Qv{pi,ni) 

_ Tr [Q(-mj ±Pihj +PiT mjhi)(l =F 7s)g(^/ TPfhf +Pf± m f n f )] ( 34 ) 
v(pf,nf)Qu(pi,ni) 

_ Tr [Qjrrij ±pjhi +pj± m^)(l =f 7s)g(~^/ TjP/jy +_P/ =F wy^/)j ( 35 ) 
8 \A(?Pi) ± nkfani)] [(QPf) T "v(?™/)] 
v(pf,nf)Qv(pi,m) 

_ Tr [Q(-mi ± pihi+pi =f m i n i )(l =F 7s)<?(-"V TPffif+p f =F m/n/)] ( 36 ) 



8y[(SPi) =F m *(^)] [(«P/) =F m /(? n /)] 

As noted above, Q is the matrix operator which characterizes the interaction. It is a 
linear combination of the products of the Dirac 7-matrices (or their contractions with 
4- vectors). 



If in the course of numerical calculations the denominators in (|26|) - ( p6|) become 
zero for some values of vectors pi, n« and pf, nf, then it suffice to change the values of 
arbitrary 4-vectors q or r, I contained in these formulae (simultaneously for all the lines 
of the diagrams being considered). 

To simplify the calculations one can substitute the vectors of the problem instead of 

arbitrary vectors q or r, I but this may cause an ambiguity of the type - which can not 
be resolved. 

The listed formulae are simple enough and may be used for creating the computer pro- 
grams for automatic calculation of the amplitudes. Besides the presence of the multipliers 
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(1 ±75) in all the formulae allows us to use the formulae of the Fierz transformations (see, 
e.g. [pip, if it is necessary, to simplify the calculations: 



[(1 ± 75)74. [(1 =F 75)7% = 2[1 ± 75]a[l =F 7fe]« (37) 

[(1 ± 75)7,4 [(1 ± 7fe)y]« = " [(1 ± TshAu K 1 ± TbJ^ljy (38) 
(i, j, k, I are indices that label the components of 4 x 4-matrices). 

6 Methods of the calculation of the amplitudes based 
on multiplication of 7-matrices and bispinors which 
are written in the matrix form and methods based 
on transformation of bispinors 

Throughout this Section we use the chiral representation of the 7-matrices: 









' a k ) 




(? 


0) 


,7*=( 


-a k j 


,75 



-J O 

O I 



s--{\\y--{°i)^--{\\y--{\\)'°--{ii)- 

In this representation the projection operator for a massive particle has the form: 
"Pip.n) = i(p + m)(l + 7 5 n) = -• 

/ m + [pn] 0z — i\pn] xy — p^n_ p_ — mn_ — + ron^ \ 

p + n± — p±n + m — [pn] 0z + i[pn] xy —p± + mn± p + — mn + 

p + + mn + p* ± + mn* ± m + [pn]o z + i[pn] xy p + ri]_ — p* ± n + 

V Ri_ + mn± p_ + mn_ — p±n_ m — [pn] 0z — i[pn] xy j 

(39) 

where for any 4-vectors a and b 

a± = a ± a z , a ± = a.,, + m y , [a6] 0z = a b z - a z b , [ab] xy = a x b y - a y b x , 



p = E = ym 2 + p 
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The projection operator for a massive antiparticle: 

1 1 
'Pa{p,n) = -(p-m)(l + j 5 h) = -• 

/ — m + [pn}Q Z — i[pn] xy p_n\ — p* L n_ p_ + mn_ —p± — rnn* ± \ 

p + n± — p±n + — m — [pn] 0z + i[pn] xy —p^_ — mn± p + + mn + 

p + — mn + p* ± — mn* ± —m + [pn] 0z + i[pn] xy p+n* ± — p*j_n + 

\ p± — mn± p_ — mn_ p~n± — p±n_ —m — [pn]o z — i[pn] xy j 

(40) 

As noted above, the calculations of the amplitudes by multiplication of 7-matrices and 
bispinors which are written in the matrix form are reduced to the general scheme @ . Let 
us illustrate this by several examples. 



1.1. Let us consider 



u(p, n) 



V(p,n) 



Tr[V(p,n)V(r,l)} 



1 



f 1 \ 



1 

V J 



2^(Po + m){\ +n z ) -p z n 
where = (1,0,0,0) , P = (0,0, 0,1) 

P(p,n) 



f P- + m(l - ra_) + [pn] 0z - i[pn} xy \ 
(m + p+)n ± — p ± (l + n + ) 
p + + m(l + n+) + [pn} 0z + i[pn] xy 
\ [m + p-)n± + p±(l — nJ) 



u'(p, n) 



Tr [V{p, n)P{r, V) 



1 


V 1 / 



2\/(Po +m)(l -n z ) +p z n 



( (m + p„)n^_ — p* ± (l + nJ) 
p + + m(l - n + ) - [pn] 0z + i[pn} xy 
(m + p + )n* ± +p* L (l- n + ) 
V p- + m(l + n_) - [pn] 02 - ifem]^ / 



where = (1,0,0,0) , I" 1 = (0, 0, 0, -1) . Note that 

u'(p,n) = e lifi(p ' n) u{p,n) , 

where 

A 9 (p,n) _ (go + m ) n l ~ Pl n o 



(po + m)(l -n*) +p z n J(p + m)(l + ra*) -p 2 n 



(41) 



(42) 



(43) 



(44) 
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Similarly: 



v(p, n) 



Tr[V a {p,n)V{r,l)} 



1 



2 J (p - m)(l + n z ) - p z n 



( 1 \ 


1 

V o ) 



( P-- m(l - n_) + [pn] 0z - i[pn] xy ^ 
(— m + p+)ri_L - + n + ) 
p + - m(l + + [pn] 0z + i fern] a,,, 
V (— m + p~)n± +px(l — n ) 



v'(p, n) 



V a (p,n) 



Tr [V a {p,n)V{r, I')} 



1 


V i J 



i 



2^(p - m)(l - n z ) +p 2 n 
Note that 



/ (— m + — + n_) \ 

p + - m(l - n + ) - [pn} 0z + i[pn} xy 
(-m + p + )n* ± +p* ± (l - n + ) 
V p- - 777,(1 + - [pn] 0z - ifpn]^ / 



where 



,i<j)(p,n) 



v'(p,n) = e i(f>(p ' n) v(p,n) , 
(p -m)n* ± - p* ± n 



(p - m)(l - n z ) +p z n J(p - m)(l + n*) - p 2 n 



(45) 



(46) 



(47) 



(48 



Obviously the bispinors ([II]) - (fl2p for particles and the bispinors (fI5p - ([IB) for 
antiparticles satisfy as the Dirac equation so also the equation for the vector n 
determining the axis of the spin projections, since when constructing the bispinors 
there were used the corresponding projection operators (|39|), fl40|). Besides one can 
easily verify that these bispinors satisfy the normalization conditions 

u(p, n)u(p, n) = u(p, n)u'(p, n) = 2m , 



v(p, n)v(p, n) = v'(p, n)v'(p, n) = —2m . 
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It follows from (gl]): 

u{p f , n f )Qu(pi, Hi) = Tr [Qu{p h n^uipf, n f )\ 



Tr 



Q 



V(pi,ni 



( 1 1 ^ 



10 10 

V o o o o J 



V{p f ,n f ) 



Tr[V(p f ,n f )V(r, I)] 



!Tr[P(pi,m)P(r,l)] 

Tr [QP( Pi ,n t )P(r, l)P(p f ,n f )] 
Tr [P( Pi , m)P{r, I)] Tr [P{p f , n f )P(r, /)] 

i.e. the calculation of the amplitude in this case reduces to the general scheme 
if 

Z = P(r,/) = i(l + 7 °)(l-757 3 ) • 

Similarly 

u'(pf,nf)Qu'(pi,ni) 



(49) 



Tr 



Q 



P(Pi,ni) 



^0000^ 
10 1 


V o i o i J 



n, 



Here 



Tr[P(p h n t )P(r, I')] 

Tr [QP(Pi,ni)P(r, l')P(p f ,n f )] 
Tr [Pfa, m)P(r, V)] Tr [P(p f , n f )P(r, V)] 

Z = P(r,r) = i(l + 7 °)(l + 757 3 ) 



Tr[P(p f ,n f )P(r, I') 



(50) 



And, as it follows from (f|^ 
differ by the phase factor. 

Similarly, for the amplitude 

we have 
and so on. 

Note that until the present time the bispinors of the general form fl4ip - ( fl"2|) , (|45D - 
( [46|) have not been used for the calculation of the amplitudes of processes, however 
their special forms described below are widely used. 



- (0), the expressions obtained for the amplitudes 

u(p f ,rif)Qu'(pi,ni) 
Z = P(r,l)-e i(p(Pi ' ni) 
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1.2. Let us consider the form of the bispinors ([IT]) - (f^2"D, (^Sj) - ( p)| ) in the situation 
when the polarization state of particles and antiparticles is the helicity, i.e. 



A 



Po 



In this case the bispinors take the form 
u(p, +) 



2^(p + m) \p\ (\p\+p z 



A = ±1 



f (po + m- \p\){\p\ +p z ) \ 
(p + m-\p\) p ± 
(p + m+ \p\) (\p\ +p z ) 
V (p + m + \p\) p ± 



( y jpo ~ \p\ (\P\ +Pz) ^ 
Po ~ \p\ PL 



2|p1 +p z 



PO + \P\ (\P\ +Pz 



V ypo + \p\ p± 



J 



In (|51p we used the identity: 



Po + m ± \p\ 
2(p + m) 



Further 



u(p, -) 



2^J(p + m) \p\ (\p\ -p z 



( (p + m+ \p\) (\p\ -p z ) \ 
(po + m+\p\) {-p±) 
(p + m- \p\) (\p\ -p z ) 
V (p + m - \p\) (-p ± ) j 



' \ jpo + \p\ (IpI -p*) ^ 
po + IpI (-p±) 



V 2 IpKIpI-a 



Po - |p| {\P\ ~Pz) 

\ \jpo - \P\ {-Pa) ) 



(51) 



(52) 
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u'(p, +) 



2j(Po + m) \p\ (\p\ -p z ) 



( (po + m- \p\) p* ± ^ 
(po + m- \p\)(\p\ -p z ) 
(p + m+\p\) p* ± 
V (p + m+ \p\)(\p\ -p z ) j 



1 



2\p\ (\p\ -p z 



I \ /po - IpI Pi \ 

Po ~ \P\ (IpI ~Pz) 

Po + IpI Pi 

V \JP0+\P\ (IpI - Pz) J 



u(p, 



2\/{Po + m) \p\ (\p\ +p z 



( (p + m+ \p\) (-p* ± ) \ 
(p + m+ \p\) (\p\ +p z ) 
(p + m- \p\) (-p* ± ) 
V (po + m- \p\) (\p\ +p z ) J 



2\p\(\p\+ Pz 



( \ Jpq + IpI (~p!) ^ 
Po + IpI (IpI +Pz) 

Po - \P\ (-Pi) 
V >Jpo ~ \p\ (IpI +Pz) J 



As this takes place we have 

u'(p, ±) = ±u(p, ±) 



PI 



pi + pi 



±u{p, ±) 



pi 
Ipj-I 



Further 



v{p, +) 



2 V (Po - m) \p\ (\p\ - p 
( 



( (p -m+ \p\) (\p\ -p z ) \ 
{pa-m+\p\) (-p±) 
(p -m - \p\)(\p\ -p z ) 
V (po-m- \p\) (-p ± ) J 



2 \p\ (\p\ -p z 



Po+ IpI (IpI \ 

\Jp o + \p\ {-Pa.) 

-\ /pq - IpI (IpI -p*) 
V - \/po - IpI (-p±) / 



In 



we used the identity: 

Po - m ± |p1 
2(p - m) 



±vpo ± IpI 
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v(p,-) 



2\J{Po ~ m) \p\ (\p\ + p z ) 



( (po-m - \p\)(\p\ +Vz) \ 
(p -m-\p\) p ± 
(p -m+ \p\) (\p\ +p z ) 
V (p -m + \p\) p± I 



1 



2|p1 {\p\ +p z 



f -y /po - \p\ (\P\+Pz) N 

- ypo - \ p\ p± 

\ Jpo + \P\ (\P\ +Pz 
\ \Jpo + \P\ P± 



2^{p -m) \p\ (\p\ +p 

( ■ Ipo + \P\ (-Pi) 

P o + |p| (\ p\ + Pz 
Po - \p\ {-pD 



J 

( (p -m + \p\) (-p* ± ) \ 
(p -m + \p\) (\p\ +p z ) 
(po-m - \p\) (-p* ± ) 
V (po-m - \p\) (\p\ +p z ) j 



2\p\(\p\+p z 



\ 



V -\/P0 - \P\ (\P\ +Pz) ) 



v'(p, 



2y/{po-m) \p\ (\p\ -p z 



\ 



( (p -m - \p\) p* ± 
(p -m- \p\)(\p\-p 
(po-m+\p\) p* ± 
V (po -m + |p|) -p z ) J 



2\p\ (\p\ -p z 



( -\ Jpo ~ \P\ Pi ^ 
\fpo ~ \p \ (\p\ - Pz) 
Po + \p\ Pi 



V \JP0+\P\ (|P1 - Pz) J 



As this takes place we obtain 

v'fa ±) = Tv(p, ±) 



PI 



Pl+P 2 y 



Pi 
\P±\ 



(56) 



(57) 



(5f 



(The signs ± in notation of the bispinors correspond to the helicity of particles or 
antiparticles.) 

The bispinors (|5T| ) - fl5"8|) have been considered in a number of papers. 
In particular, in p6f the bispinors (|51~1) , (|52|) were considered. 
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In [27 1 the bispinors flSTp, Q5l|) were considered. 

In [^8] it is proposed to use the bispinors ([5TD, fl5"4|), (|56|), (|57|) [in this paper the 
bispinors (|56|) , (|57|) have the opposite sign]. 

In [^9| it is proposed to use the bispinors (|5lD - (|58D [in this paper the bispinors 
([55]), ( |5TD have the opposite sign]. 

In |3(J the bispinors (|51|), (p4[), (|56"D, (|57|) are used [in this paper the bispinor (|57|) 
has the opposite sign]. 

In |TI| there are presented the expressions for the amplitudes that can be obtained 
by means of the formulae (|52"D, (|53|), (|56|), (|57|) written in the center-of-mass frame 
of the fermions considered. 

The bispinors described in this subsection are also widely used in calculations others 
authors. 



1.3. Let us consider a special form of the bispinors (|5T 

u+ (p) = v_(p) -■ 



( \ 




P+ 
\P± J 



|)atm = (|pl = p ) : 

(59) 



u_(p) = v+(p) 



u' + (p) = v'_(p) 



fpZ 



u'_{p) = v'Jp) 



'Pa 



( P- \ 
-P± 


V o j 
( \ 



Pi 

V p- ) 

( -Pi \ 

p+ 



V o j 



(60) 



(61) 



(62) 



In a number of papers it has been proposed to use the bispinors ( [59D - (162|). 
In |^| it was proposed to use the bispinors ([5T?D, (|B2]). 

The method proposed in |33|] is reduced to the application of the bispinors of the 
form (|60|), (|6l|) [bispinor (|60|) was used with the additional factor i]. 

In [^4 it was proposed to use the bispinors of the form (^), (|62] ) [here the bispinor 
(|62"D was used with the opposite sign]. 
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The method proposed in [f55|] reduces to the application of a special form of the 
bispinors (|59|) - (^2]) at p y = . 

These methods are widely used in calculations. 



1.4. Foif] 



„ A p x p z p,,p z p : 
n ^ = — [p z , : , — ,m + 



m po + m po + m 
the formulae (0), (g2|), (g^), (g§ imply 

1 

u(p, +) = 



Po + m' 



A = ±1 



2(p + m) 



u'(p, 



v(p, 



v'(p, +) 



1 



2(p + m) 



2(p Q + m) 



1 

2(p + m) 



/ p_ + m \ 

-P± 
p + + m 

V p± ) 

( -rf. \ 

p + + m 
p* ± 

V p_ + m J 

( — p_ — m \ 

P± 
p + + m 

V P± / 

p + + m 

-Pi 

V — p~ + m J 



(63) 



(64) 



(65) 



(66) 



The bispinors ( |63"D - ( |66j ) are considered in [Q, j37|, (in |37j and J^] the 
bispinor (|66| ) has the opposite sign). 

In [|39| the bispinors (|63|), ( |54]) are used for a special case when the particles are 
moving on xy plane, i.e. p M = (po,p x ,P y , 0) and n^ 1 = A(0,0,0, 1) , A = ±1 . 

1.5. Foif| 



A ?n 

™ M = — (Po : 

m Po+Pz 



Px,Py,Pz + 



Po+Pz 



, A = ±1 



the formulae (0), ((H), ©, © imply 



u(p, +) 



/ m \ 


P+ 

v p± y 



(67) 



3 In the rest frame of a particle (antiparticle) the vector takes the form A(0,0,0, 1) 

4 In the rest frame of a particle (antiparticle) the vector n M takes the form A(0, 0, 0, 1) 
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u'(p, -) 



rp-\ 



v{p, 



fp-\ 



\ 



v'{p, +) 



( -Pi \ 

P^ 


V m J 

( — m 


P+ 

V p± J 

f-pl\ 

o 

V — m ) 



(6* 



(69) 



(70) 



The bispinors (§7|) - 
bispinors take the form 



are used in [40 , [521. Note that in the massless limit these 



2. In the papers f4lfl , f42| , |43|] it is proposed to use the following bispinors for massive 
Dirac particles (antiparticles), the polarization state of which is the helicity: 



«"(p, +) 



2 \p\ | PJ 



Po 



- |p1v(ipi + p 2 )p! \ 



PO - \p\J(\p\ ~Pz)P± 



Po + \p\ y (\p\ +Pz)p*± 

\ \JP0+\P\\J(\P\ ~Pz)P± ) 



(71) 



/ _ 



u"( P , 



2 \p\ | Pj 



+ Ip1\/(Ip1 -p*)p1 \ 



\ /g0 + \P\ \ J(\P\ +Pz)P~ 
-\J P0-\P\ ^ J{\P\ -Pz)P*± 

yjpo ~ \P\>J(\p\ +Pz)p± ) 



(72) 



v"(p,+) 



2 \p\ | PJ 



+ IpV(Ip1 -p*)p1 \ 



po + IpVOpI + p^)p± 



-Vpo - IpIvGpI 



(73) 



</'(p, -) 



2 |pl | Pj 



po - \p\y(\p\ +p z )p± J 

1 \ /po - \p\ \ fi\p\ +Pz)p*± N 

\ /p0 ~ Ipj \ /(|P1 ~Pz)p± 

y/ po+ \p\ \ /(\p\ +Pz)pj 

Po + IpI v^IpI -Pz)p±_ J 



(74) 
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(In [[TJ] the bispinors for antiparticles have the opposite sign, and in [53] the bispinors 
for antiparticles have an additional factor i). 

As this takes place we have [see (|5l|) - (^), ( 



u"(p,±) = ±u(p,±) 



Pi 



v"(p,±) = ±v(p, ±) 



pi 
\p±\ 



Thus if we use the bispinors in form ( |TTD - ( |T4] ) then the calculation of the amplitudes 
of processes is reduced to the general scheme (0). For example for the amplitude 



u"(p f ,+)Qu"(pi,+) 



we have 



V(r,l) 



(Pi)i 



(Pfl 

\\(Pi)±\ \\(Pfl 



where = (1,0,0,0) , P = (0,0, 0,1) . 
3. In the paper [Q it is proposed to use the bispinor for a massive particle in the form: 



u(p, n) 



V(p, n) 



J (pq) — m{nq) 



( \ 


Q+ 

\q± J 



where 

— {Qo = 1 + |a| 2 , q x = 2Rea, q y = 2Ima, q z = 1 — |a| 2 ) is a massless 4-vector, 
a is an arbitrary complex number and V(p,n) has the form (|39|). 

In this paper it was shown that 

Tr [QP(pi,7H)V + (q)V(p f ,n f )] 



u(p f ,rif)Qu(pi,ni 



Tr [V( Pi , ni)V + {q)] Tr [V(p f , n f )V+(q)} 



where V+(q) = -(1 + 7 5 )<? . 

So the methods of the calculation of amplitudes discussed above in this Section are 
described by the general scheme ([|). Moreover, most of them are the special cases of 
the covariant method proposed in Section 4. However the defects of these methods are 
complicated calculations, bulky and noncovariant form of their results. 

In conclusion let us consider two methods based on bispinor transformation. These 
methods are similar to the ones discussed above, but do not use the explicit form of the 
bispinors and projection operators. 
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4. In the paper it is proposed to transform the bispinor for a massive particle as 
follows 



u(p, n) 



V(p, n) 



(pq) + m(nq) 



u-(q) 



V(p, n) 



Tr[P(p,n)V-(q)} 



(75) 



In fact, the expression for the bispinor is given by ( [75] ) up to a phase factor. Really, 
the right-hand side of Eq. (|75D can be written as 



V(p,n) 



Tr[V(p, n)V-{q)} 



u(p,n)u(p,n) u(p,n)u_(q) 
u -(q) — T~, n — rr\ u -(Q) = ~r~t \ — 7~V\ u \Pi n ) ■ 



\u(p,n)u-(q)\ 



\u(p, n)u-(q)\ 



Further 



(p + m)(l + 7 5 n) (p + m)q(p + m)(l + 7 5 n) 
u(p,n) = u -(q) = . — ""-(<?) 



2 V (pq) + m(nq) 



4(p?)v (PQ) + m(nq) 



(p + m) |(1 - 75)5 + Z|(l - 75)^/ P(p, n 



-it. 



4(pq)y(pq) + m(nq) 
(p + m) \V-(q) + lV-(q)l\ V(p, n 



(Q) 



■ u -(q) 

^(vq)\J(pq) + m(nq) 
(p + m) {Tr [P-(g)P(p, n)] + /Tr [p_(g)fP(p, n)] } 



(76) 



M -0?) 
^(pq)\I(pq) + "^(rag) 

(p + m) | [(pq) + m(nq)} + [(pi) (nq) - (pg) (nZ) - is(p, n, q, I)] l\ 



^(pq)\/(pq) + m(nq) 



u-(q) , 



where q is an arbitrary massless 4- vector; I is an arbitrary 4- vector such that 
I 2 = — 1 , ql — [to obtain fl76l), we used a variation of the identity (|21~D 

P±( 9 )Ai±(g) = rr[P±(g)A] ■ u±(g) 

where A is an arbitrary matrix operator]. 

As it follows from (75), the above presented transformation of the bispinor leads to 
the method of the calculation of the amplitudes given by a special case of formula 
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5. In the paper it is proposed to transform the bispinor for a massive particle as 
follows 



. , u + (q p )u_(q m ) , 
u(p,n) = + ^ ' W J u+((f) + U-(q m ) , 



m 



(77) 



where 



1 1 m 2 

q p = -(p + mn) , if" = -(p - mn) , (q p ) 2 = (q m ? = , (q p q m ) = — ■ 

This transformation can be obtained by the method similar to the one discussed in 
previous subsection. Really, 



u(p, n) 



V{p, n) 



Tr[V{p,n)V-{q m )} 



u-{q m ) 



q p + m i(l +7 5 )g p + m 
u-(q ) = 



m 



m 



u-{q r 



u+(q p )u + {q p ) 



+ 1 



in 



u + {q p )u.{q m ) 



u + (q p )+u4q r > 



Because of this, the transformation of the bispinor considered leads to the method 
of the calculation of the amplitudes given by a special case of formula (j33|). 

However, to calculate u + (q p )u^(q m ) the authors of [^6| use the method which was 
proposed in the papers [El - |ES|. (This method was discussed in subsection 10 of 



Section 2.) In this case we have up to a phase factor [see ([28|)j 

Tr[(l - ^)q m lqq p ] 



u + (q p )u4q m ) 



4j(qq p )(qq m ) 



W)y + i{q P U 



(q m )o - (<f 



--[(Ov + iGTU 



(q p )o - (q p )x 

\ (q m ) ~ (q m ) x 



l \ (q p )o-(q p ), 

[where P = (0, 0, 1, 0) , q» = (1, 1, 0, 0) ]. 

As a result, the first term in ([F?]) obtains the phase factor [see (10)] 

U-{q m )lqu + (q p ) 
\u_(q m )lqu + (q p )\ 

and the bispinor under consideration does not satisfy the Dirac equation. 

Note that a special case of the transformation considered, namely the situation when 
a 4-vector, specifying the axis of the spin projections, has form 

1 m 



n = A 



— V 

m 



(pq Tl 



: ( 1 



A = ±1 



was proposed in EH . 
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6. Incorrect methods of calculation of the amplitudes 

In JT7 for deriving of expressions for amplitudes incorrect identity (15), writting with- 
out any substantiation, is used 



cc? = —(ss') + 



jps')(p's) 
(j>p') + mm' 



(15) 



where s = — , s = c H ■ p ; s = — - , s = c + — -p 

m po + m m p + m 



Authors of |48| propose to use bispinors for massive Dirac particles in the form 



u+(p) 



im 



2(pg)g_ V ?± 
( kl \ 



u-ip) 



m 



2(pq)q_ \ V- 



where 



p = k + 



m 



z 2 2 n 

g, p = m , g = U , 



2(pq) 



m 



k = p- rq ; A; 2 = 0. 
It is stated, that g is an arbitrary massless vector. However from elementary requirement 

u+(ph 5 u+(p) = u_(p)75«-(p) = 

we have 

q x Px 



Qo -q z Po- Pz 

i.e. vector q is not arbitrary. Besides, 4-vector, which determines the axes of the spin 
projections, is not defined in this paper. 



In paper JEJ there are mistakes in phase factors of coefficients Cw, given by equations 
10) therein. As a consequence method gives incorrect results. For example, using this 
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method we have 

" + 



Hp, +)jsu(p', +)] [u(p, +)jsu(p', +)}* 



\{pp' + mm') 2 J pi + PyP 2 

2 



{(POPO - \p\W\ + m m')\jpl+P 2 y [Po{0) ~ P'of 



+ (PoPo + \p[W\ + mm ')^Pl+P 2 y[\p\W\ ~ W)} 2 m 2 
-2(mp' + m'p )p x p G {0) - p' f [\p\W\ - [0)\ m} 
However, calculation by the classical method gives another result: 

|«(p,+)75ii(p',+)p 



2 [mm 1 - pop'o + \p\\f\][ffl) + 



7 The computer programs for the calculation of the 
amplitudes 

A number of computer programs for the calculation of the amplitudes of processes with 
polarized Dirac particles (antiparticles) has been developed in the last few years. However 
all of them have these or those disadvantages and limitations. In particular, they do 
not allow us to calculate the amplitudes of processes involving massive Dirac particles 
(antiparticles) with arbitrary polarization states. Let us consider these programs briefly. 



1. The program CHANEL (see |5(|) is a part of the computer system GRACE (see [pl| ). 
This program allows us to calculate the amplitudes of processes with massive Dirac 
particles (antiparticles), the polarization state of which is the helicity, i.e. the 4- 
vectors specifying the axes of the spin projection are 

[where p^ = (po,p) is the 4-momentum of the corresponding particle (antiparticle)]. 

The program uses the method of calculation, which is analogous to the one consid- 
ered in subsection 5 of Section 6: 

(a) 

A = +1 
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u(p,A = +l) 



V(p,X = +l) 



Tr[V(p,\ = +l)V + (q P + i)} 



m 



|(1 - 75)9+1 + m 



(b) 



u-te)«+(£i) 



W _(9+\)+m + (9+i) 



(78 



2b1 



A = -1 
(|P|,-P) , Gft)' 



u{p,\= -1) 



P(p,A = -l) 



Tr[P(p,A = -l)P-(^ m i)] 



2b1 

«-(9 m l) = 



9_1 + 



77? 



m 



«-(<&) 



777 



777 



«+(£!)+ 



Introducing the notation 



n — rP — n m n — n m — rP 

Qi — 9+1 — Q-i > 92 — q+i — q~i 



and combining ( [75[ ) and ( 75 ) we obtain 



u(p, A) = M_ A (<? 2 ) + «A • 



777 



Further we have [up to phase factor - see also fl28 



w-a(?2)wa(9i) - A- 
Ap y - 



qiqiq 2 (l + A7 5 ) 



4\/ (9?2) 



Am- 



Tr 



g/f/ p (l - A75) 



4 v(?r) 2 - (?y< 



Py 2 + p. 2 



Here 



(79) 



30) 
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r» = (1,0,0,0), //=(0,-^), 

IP] 

q, I are 4-vectors such that 

q 2 = 0, l 2 = -l, (ql) = 

[in this case one chooses q^ 1 = (1, 1,0,0) , F = (0,0, 1,0)]. 
However in this case the first term in (|80"D obtains the phase factor 



\u x (qi)qlu-x(q 2 )\ 

As a result, the bispinor (ROT) does not satisfy the Dirac equation. 



The program COMPUTE (see |52|) use the method of the calculation of the ampli- 
tudes of processes with massive Dirac particles (antiparticles) proposed in |40| and 
considered in subsection 1.5 of Section 6. 

Remind that in this case the polarization state of particles (antiparticles) is defined 
by the vectors 



n 



A , m m . 

— [Po ; ,Px,P y ,Pz H ■ ) , A = ±1 

m Po+Pz Po+ Pz 



[where = {po,Px,Py,Pz) is the 4-momentum of the corresponding particle (an- 
tiparticle)]. 

3. The program HELAS (see ||53|| ) is a part of the computer system MadGraph (see [p4"|). 

It calculates the amplitudes of processes with massive Dirac particles (antiparticles), 
the polarization state of which is the helicity by the method proposed in |26|| . This 
method is considered in subsection 1.2 of Section 6. 

4. There are two variants for the program HIP. The first variant (see [j5"5|] ) calculates the 

amplitudes of processes with massless Dirac particles (antiparticles) by the method 
proposed in []2S| and considered in subsection 10 of Section 2. 

The second one (see fl56|) calculates the amplitudes of processes with massive Dirac 
particles (antiparticles), the polarization state of which is the helicity, by the method 
proposed in |26[ and considered in subsection 1.2 of Section 6. 



5. The program SPINORP (see f57j) calculates the amplitudes of processes with massless 
Dirac particles (antiparticles) by the method proposed in [Q and considered in 
subsection 1.3 of Section 6. 
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6. The program WPHAST 1.0 (see J5BH ) calculates the amplitudes of processes involving 



massive Dirac particles (antiparticles) by the method proposed in |16| and considered 
in subsections 5.2 and 10 of Section 2. 

Recall that in this case the polarization state of particles (antiparticles) is defined 
by the vectors 

1 m 



n = A 



— p- 7-^9 

m ypq) 



A = ±1 



[where p is the 4-momentum of the corresponding particle (antiparticle) and q is 
an arbitrary 4-vector such that q 2 = and in numerical calculations one chooses 
<f = (1,1,0,0)]. 

7. In there was proposed the program for the calculation of the amplitudes of pro- 



cesses with massive Dirac particles (antiparticles) by the method proposed in [KH] 
and considered in subsection 1.4 of Section 6. 

Recall that in this case the polarization state of particles (antiparticles) is defined 
by the vectors 

„„ = A fc ,J?* J^, m+ _d_), A=±1 
m p + m p + m p + m 

[where p^ = (po,Px,P y ,Pz) is the 4-momentum of the corresponding particle (an- 
tiparticle)]. 

However, in addition to programs described above some authors use the computer 
algebra systems to calculate the amplitudes of processes for concrete problems. For ex- 



ample, the author of |28[ uses the system MAPLE for the calculation of the amplitudes of 
processes involving massive Dirac particles (antiparticles), the polarization state of which 
is the helicity, by the method considered in subsection 1.2 of Section 6. 

Thus, at present there are no computer programs for the calculation of the ampli- 
tudes of processes involving massive Dirac particles (antiparticles) which have arbitrary 
polarization states. The programs of this sort may be created with help of the formulae 
presented in Section 5. 
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